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Abstract 

We investigate on the lattice the Yukawa models in 2 dimensions with Z(2) and U(l) 
symmetries. These models reduce to the usual and chiral Gross-Neveu models, respectively, 
when the kinetic and the selfcoupling terms of the scalar field are turned off. The numerical 
data and mean field arguments suggest that, at least for some range of the scalar field hopping 
parameter, fermion mass is dynamically generated for arbitrarily weak Yukawa coupling. The 
models are asymptotically free in this coupling, like the Gross-Neveu models, even when the 
scalar quartic selfcoupling is strong. 
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1 Yukawa models in 2 dimensions 



The 2d Yukawa models (Y2) have received little attention since the rigorous establishment 
of some of their fundamental field theoretical properties in the seventies (for a summary see 
ref. The superrenormalizability of the Yukawa coupling and the arbitrariness of the 

selfcoupling terms of the dimensionless scalar field in 2d might have led to a feeling that the 
Y2 models have little relevance for the 4d field theories. This is to be compared with the 
continuing interest in the 2d Gross-Neveu (GN) and nonlinear a (NLc) model investigations, 
shown e.g. by the recent exact mass gap calculations in these models @, ||. 

However, the Y2 models can be chosen so that both the GN and NLcr models are their 
special cases. This is most obvious on the lattice. For example, the action of the Z(2) 
symmetric Y2 model on the lattice can be chosen in the form 

s = -2«x)^«+M+E^+^E(^- 1 ) a 

X,fl X X 

+ y:c#:+i/E&c- (1.1) 

x,a x,a 

Here we have introduced N "naive" Dirac fermion fields ip a , a = 1,...,N, on the lattice which, 
due to the fermion doubling, describe Np = 4N Dirac fermions of zero bare mass. All the 
fields and the couplings are made dimensionless by the appropriate rescaling with the lattice 
constant a, x enumerates lattice sites and is the lattice derivative. The action for the 
U(l) symmetric Y2 model which we study as well has the form similar to fll.lj ) with two- 
component field (j) x = ((/)]., (j) 2 x ) instead of cj) x and the Yukawa term y Y, x ,a ^"(^l + i^tlp)^ 
(where 7p is the 2d analog of 75). 

At k = A = the action fll.lf ) describes the Z(2)-symmetric GN model in the auxiliary 
scalar field representation of the 4-fermion coupling. The Yukawa coupling y is related to 
the usual GN coupling g by y = V2g. On the other hand, at A = 00 and y = the action 
( |1 . 1| ) describes the Ising model. The U(l) symmetric Y 2 model in similar cases reduces to 
the chiral GN model or to the XY model, i.e. the NLa model with U(l) symmetry, and 
similarly for other symmetry groups. The Y2 models thus interpolate between the GN and 
spin or NLa models. 

On the lattice, by choosing the above formulation of the scalar field sector using the 
hopping parameter k, the kinetic term can be turned on or off gradually, illustrating the 
smoothness of the transition from an auxiliary to a dynamical scalar field. As the numerical 
study of lattice Yukawa models in 4d showed, the physical observables behave continuously 
with k in the vicinity of k = including a region of negative k @, |5|, This fact elucidates 
the relation between the Yukawa and four-fermion theories found in the continuum in the 
leading order in 1/Np expansion ||. |7j. 

Motivated by these considerations and by the recent discussion of a relationship between 
the Nambu-Jona-Lasinio type four-fermion theories and the Standard Model f|, [7| (for a 
recent review see ref. M), we address here the question to what extent the Z(2) and U(l) 
Y 2 models still possess the most interesting and important properties of the GN models 
[[], |HJ, namely the asymptotic freedom of the Yukawa coupling y, the dynamical fermion 
mass generation and, in the case of the Z(2) model, the dynamical symmetry breaking. 

2 Expected scaling properties 

In the GN models the basic scaling properties at y — > can be derived by means of the 1/Np 
expansion. For Y 2 models this expansion is applicable only for small A , strictly speaking 



for A = 0(1/ Np) @, 0]- It gives results very similar to those for the GN models, provided 
k < k c (X). Here k c (A) is the critical line of the scalar model at y — (i.e. without fermions) 
which for — k c (A) < k < k c (A) is in the high-temperature phase. The fermion mass in lattice 
units amp is expected at fixed k and A to scale with y as 



amp oc exp 



1 a 2 m\ 1 



(2.1) 



Here /?o is the first coefficient of the /^-function of the continuum GN model with Np flavours 
of Dirac fermions, 

A = -^ (Z(2)), (3 = ^ (U(l)), (2.2) 

am^ is the 0-field mass in lattice units and its renormalization constant at y — 0. For 
A = we have at y = the free scalar field theory with 

aS < = t 1 - ^oy) \ • z * = h • < 2 - 3 > 

and k c (A = 0) = 1/4. Note that for small A the ratio Z ( j ) /a 2 m'l appearing in eq. ( p.l| ) is the 
scalar field propagator at zero momentum, i.e. the susceptibility x °f the pure scalar model 
at the same values of A and k. Although a 2 m\ and are not well defined separately at 
k < 0, the ratio x is- Thus the scaling law fl2.1|) is sensible also at k < 0. We shall see below 
that for small A the mean field approximation (MFA) gives the same results as the 1/Np 
expansion. 

For large A the 1/Np expansion is a priori not applicable. However, as we shall argue in 
the next section, the MFA can be applied to study the Y2 models at all A > 0. The effective 
interaction produced by fermions is very nonlocal and favours ferromagnetic ordering of the 
field 4> x . Such an interaction can be well described by some effective mean field H acting on 
(f) x at each site. Thus using the MFA we reduce the model to a pure scalar model with only 
local interactions and the external field H. Such a model can be easily studied e.g. by high 
temperature expansion or Monte Carlo (MC) simulation. 

The results we obtain using this approach are as follows: 

(i) At any nonnegative A including A = 00 the Y 2 model with Z(2) symmetry is in the 
phase with broken symmetry (((f)) 7^ 0) for arbitrarily small y. Similarly, the corresponding 
U(l) model is in the spin wave phase (analogous to the low temperature phase of the XY 
model). In both cases fermion mass is dynamically generated for arbitrarily small y and 
vanishes only at y — 0. This means that the Yukawa coupling is asymptotically free also for 
any nonnegative A. 

(ii) For n < k c (X) the fermion mass and, in the Z(2) case also the magnetization y(<p), 
scale with y according to 

h(n, A) 



am F oc exp 



y 2 



(2.4) 



(in) In the MFA the function h(n, A) is given by 



h{K t X) = ~ t (2.5) 
N F x 



where x is the zero-momentum scalar propagator in the pure scalar 4 model at given k 
and A. If the propagator is dominated by a pole with mass am^ and residue Z^ y then 



Thus according to the MFA the asymptotic freedom and the other mentioned properties 
of the GN models occur also for large A in the Y 2 models. Our Monte Carlo data at A = 0.5 
and A = oo support (i) and (ii) (see also ref. fll!!), whereas some deviations from (Hi) are 
observed. 



3 The mean field approximation 

Integrating over fermion variables in the partition function for the Z(2) Yukawa model ( |1 . 1| ) 
we obtain an effective scalar model with the contribution to the action from the fermion 
determinant 

S det [<P] = -Ntr\nM[<j ) \ 1 (3.1) 

where ^ 

M[4>] xz = -jpiS^+p, - 5 ZjX ^) + y<f) x 5 xz = K xz + y<pj xz . 

We know from the experience with Yukawa models in four dimensions (for a recent review 
see ref. ||) that fermions strengthen ferromagnetic ordering: the transition from disordered 
to ordered phase occurs at smaller values of k as y increases from zero. As was pointed out 



to us by E. Seiler [12] one can prove also that Sd e t[0] is minimized on the totally ordered 
configuration. 

We first give a very straightforward derivation of the announced results by the MFA 
method and then discuss its reliability. As infinitely many sites participate in the interaction 
Sdet with a given one, say x, we expect their effect on <p x to average into an effective mean 
field H whose fluctuations are negligible. So we obtain a scalar <p A model with only local 



interactions given by the first three terms in (1.1) to which the constant external field H is 



applied. The mean magnetization a of such a model is given at each k and A by the response 
function o = f(H). 

The simplest approach to calculate H is to differentiate Sdet with respect to 4> x at the 
site x and then substitute the mean value a for all 0's. This gives 

H(a) = 2Ny 2 a J ^ (E 8 ^ + (^) 2 ) • (3-2) 
Selfconsistency then requires that 

a = f(H(a)), (3.3) 

where H(a) is from ( p.2[) . For given y, k and A one can find a as a solution of this equation. 
At k < k c {X) for small a we can write 

f(H) = X H + 0(H 3 ), (3.4) 

where \ is the susceptibility of the scalar 4 model, and (|3.2|) gives 

H(a) = -Ny 2 a In — + 0(a). (3.5) 
7r ya 



Substituting ( |3.4| ) and ( |3.5| ) into ( |3.3| ) one finds that a nonzero solution exists for arbitrarily 
small y. Taking amp = yu we arrive at eqs. fl2.4|) and (|2.5|) with Np = 4iV. Only the slope x 



of the response function f(H) at the origin was needed to derive this asymptotic scaling law. 
As we shall see in Sect. 4, the eq. ( |3.3|) with the full response function can give a very good 
approximation for mp even when a is not small. 

Now let us discuss the MFA method in more detail. To get an idea of how nonlocal 
the interaction S^et is at small y and whether it can produce the effective mean field let us 
expand it in powers of y: 

S det {$\ = const + N^J2 TrK^K^fafa + 0(y 4 ), (3.6) 

xz 

where Tr is the trace over Dirac indices. The term of order y 2 produces a nonlocal interaction 
between <p x and (p z \ 

S 2 = -Ny 2 J2Jx Z <Px<Pz (3.7) 

xz 

with J xz = ^TrK~ z K~ z ( as K~ z = —K zx x \ One can find out easily that J xz > when x 
and z are separated by an odd number of links and J xz = otherwise. Thus this interaction 
is ferromagnetic. 

The MFA works well when fluctuations of the mean field at site x, H x = 2Ny 2 J2 Z Jxz<Pz, 
are small compared to (H x ). In 2 dimensions J xz falls off with \x — z\ so slowly that it 
produces the infrared logarithmic divergence 

£ J xz = l -TrK~ 2 =/^(E sin2 *v) • ( 3 - 8 ) 

or, in other words, very large number of sites contributes effectively to H x . This means that 
the mean value of H x grows with the size L of the system as Ny 2 (<p) \nL if there is some 
nonzero (4>). As one can check the fluctuations of H x around its mean value are finite at 
L — > oo and thus at large L the MFA is justified if S 2 was the only interaction. Moreover, 
the divergence in (H x ) means that such a system would be ordered at any y. This can be 
also understood if one considers the free energy of the long wavelength modes of 0. It is 
negative and diverges as InL. 

The interaction Sdet contains also other nonlocal terms which are of higher order than y 2 . 
They obviously contribute to H in the expression ( |3.2| ): the mass term amp = yo can be 
interpreted as their effect. However, one can easily realize that the interaction producing the 
mean field H{a) in ( [3. 2D almost coincides with S 2 for the distances smaller than 0(l/mp) 
and falls off exponentially at larger distances. This means that when mp becomes smaller 
the MFA works better as more sites participate effectively in the interaction. The size I/tuf 
plays the role of the infrared cutoff which controls the logarithmic divergence in eq. ( |3.5|) . 

It is also very useful to look at what happens from the point of view of the free energy 
of the long wavelength modes, i.e. effective potential. The effective potential V e ff(a) at the 
origin behaves like V e ff(a) ~ y 2 a 2 lnya (this is clear from the fact that H(a) in ( |3.2|) is 
minus the derivative of the fermion contribution to the V e ff(a)) and thus a = is a local 



maximum. Approach based on the effective potential was applied in a recent paper |n| and 
gives the same results as the MFA method, as one would expect. We think, however, that 
the MFA approach allows us to understand better the nature and the effect of the infrared 
singularity, the mechanism of the fermion mass generation in the U(l) case (see below) and 
the subtleties of both methods. 

Although our MFA method as well as the effective potential approach [y| and their 
results might look very consistent one should be cautious and not overestimate their relia- 
bility. The expression ( |3.2|) has a form of a contribution of one loop with massive fermions. 



Higher loop contributions may turn out to be important. For example, for the GN model 
(k — A = 0) the selfconsistency equation (|3~3| ) coincides with the gap equation obtained 
in the leading order in 1/Np expansion. However, we know that in the Z(2) GN case the 
coefficient h of 1/y 2 in the scaling law given by ( |2.5| ) receives 1/Np corrections (compare 
with eq. (|2.2|)). The same result for h as eq. (|2.5|) is obtained in [DJ and must receive finite 



corrections for this reason as well. The terms in V e ff of higher order in y neglected in |T3 



contain in fact powers of y 2 \n(l/ya) and are not negligible as y — > 0. Nevertheless, one can 
expect that H still has the singularity similar to ( |3.5| ) at m F — > as a consequence of the 
strong nonlocality of the interaction Sdet in 2 dimensions. One can also expect that m F ^ 
if cr 7^ 0. Then the conclusion (i) that the fermion mass is generated for arbitrarily small y 
still holds. It is natural to expect that these features do not depend on the structure of the 
local scalar interactions in ( p. . 1|) , in particular on the value of A. The predictions (it) and in 
particular (Hi) are less reliable, however. 

Let us now see how one can apply the MFA in the Y 2 model with the U(l) symmetry, 
when the symmetry cannot be broken spontaneously. It is known that this fact does not 
prevent fermions from acquiring a mass [|l(J. To understand how this happens we consider 
the interaction of the scalar field induced by the fermion with a small mass mp in one loop 
that we have already discussed. For distances between sites smaller than 0(1 /mp) one 
can neglect the mass and the interaction behaves like ( |3.7| ). Such interaction can produce 
ferromagnetic ordering on the distances smaller than 0(1/ nip) if mp is small enough. On 
the distances larger than 0(1 /mp) long wavelength fluctuations (spin waves) destroy the 
ordering as usual in two dimensions. However, in a finite volume of the linear size 0(1/ mp) 
the magnetization is nonzero and its direction is drifting slowly. 

Now let us imagine a fluctuating scalar field with such properties playing the role of 
a background field coupled to a fermion via Yukawa interaction. For small y we expect 
the fermion to acquire a mass given by amp ~ y(4>)- What ((f)) enters this formula? It is 
reasonable to expect that the fermion "feels" only the ordering on the scale of its Compton 
wavelength, so that (0) should be averaged over a volume of the size 0(1/ mp). This agrees 
with the observation made by Witten that "almost long-range order" in chiral GN models 



is sufficient to generate the fermion mass (10 



Bearing this in mind we can carry out the mean field considerations in the U(l) case 
in complete analogy with the Z(2) case. One has to realize, however, that ((f)) = a in the 
formulas is the average magnetization on the scale of 0(l/m F ). We actually see nonzero ((f)) 
in a finite volume in MC simulations which agrees with amp = y((f>) remarkably well (within 
few percent) when L « 1/amp. 

We now compare the solutions of the equations (|3.2| ) and ( |3.3| ) on finite lattices with 
numerical results and use these equations for data analysis. 



4 Numerical results 

We simulated the Z(2)-symmetric Y 2 model defined by the action ( |1.1| ) and the analogous 
model with the U(l) symmetry on the lattices L 2 = 16 2 , 32 2 and 64 2 , mostly with Np = 16. 
The expected scaling behaviour Q2.4p suggests to collect data at many ^-points at fixed values 
of A and k and to study the ^-dependence of the fermion mass. We have chosen several k 
points in the interval —0.3 < k < k c (X) with A = 0,0.5 for both Z(2) and U(l) models and 
A = oo for the U(l) modelQ. For each pair of (k, A) values and L we determined the fermion 

lr The hybrid MC algorithm for simulating dynamical fermions we are using unfortunately does not work 
in the Z(2) case at A = oo because of the discrete nature of the variables. 



mass amp at several y values such that 2/L ~ amp < 1. 



As has been reported in ref. Jl 1] in detail, the dependence of amp on y at all investigated 
(k, A) points is consistent with the expectations (i) and (it). The data analysis on finite 
lattices has been performed by means of the generalized gap equation 

h 7T ^ 1 

f = 21? £j £ sin 2 p M + (am F / S )2 ( ^ 

which gives the scaling law (|2.4j ) at L = oo and describes quite well the onset of finite size 
effects and a very slow approach to the asymptotic scaling. The values of the fit parameters 
h and s depend only on k and A but not on L. Note that eq. ( f4.1|) at s = 1 can be viewed 
as a selfconsistency equation like (|3.3|) with a linear response function whose slope is a free 



parameter. 

In the present paper we compare these results with the parameter-free predictions of the 
MFA, to see how well it can describe the data. To obtain these predictions we solve for 
given k, A and y the finite lattice version of eq. (|3.2|) and eq. ( |3.3|) numerically. On a finite 
lattice we substitute the integral in ( |3.2| ) by the discrete sum over the momenta allowed by 
the boundary conditions on the fermion fields: periodic in one and antiperiodic in another 
direction. The response function f(H) is obtained by a MC simulation of the scalar <f) A model 
in the external field. For k = the response function can be also expressed in terms of easily 
computable one dimensional integrals. So for each (n, A) pair and L we obtain a function 
amp(y) which we compare with the data. 

In the case of the Z(2) model we find only a qualitative agreement. There is a systematic 
discrepancy between the data and the MFA prediction at any n and A values. This is 
connected with the fact, discussed in the previous section, that in the scaling limit amp — > 
the one-loop formula for H(a) ( |3.2| ) does not reproduce the correct value of the /?o coefficient 
in the scaling law ( |2.1|) already for the GN model (k — 0, A = 0). We note that this happens 
in the Z(2) GN model also in the leading order of the 1/Np expansion. It would be very 
useful to include the higher loop corrections but at the moment we do not know how to do 
that at large A systematically.^ Without these corrections the simple MFA that we used 
apparently cannot describe the data in the Z(2) model sufficiently well quantitatively. 

In the chiral GN model higher loop corrections do not modify the coefficient of 1/y 2 in 
the exponent of the scaling law. Correspondingly, the agreement is much better in the U(l) 
Y 2 models. For comparison we present plots for (k, A) = (0, 0), (0.2, 0), (0, oo) and (0.2, oo). 
At A = the 1/Np expansion can be applied and in the leading order coincides with the 
MFA prediction. The data agree with this prediction well (see figs. 1 and 2) as one can 
expect if Np is large. For k = at A = 0.5 and A = oo (fig. 3) the data agree with the MFA 
prediction as well as at A = 0. However, if k ^ we observe both at A = 0.5 and A = oo a 
significant discrepancy illustrated for A = oo and k = 0.2 in fig. 4. 

We see that the MFA works not only at A = but also at large A as long as k = 
even better than one could naively expect. The discrepancy at k ^ and large A could 
be due to the crudeness of the simple MFA that we used and more systematic study is 
needed to understand it fully. In any case, as the fits with the generalized gap equation ( flTl) 



performed in ref. JTT]| demonstrate, all our data are consistent with the scaling of the general 
form Q2.4|) . Here we show such a fit in fig. 4 (dotted-dashed line). The values of the fit 
parameters in eq. (3~T) are approximately s = 0.27 and h = 0.86(7r /Npx)- We cannot draw 



firm conclusions from these numbers as they are obtained in the region where the asymptotic 



2 Somc additional complications in the subleading order in 1/Np arise in the Z(2) model even at A = 
due to the naive fermion discretization in eq. (p|). 



scaling is not yet achieved and we do not know how well the eq.( |4.1|) describes the approach 
to the asymptotic scaling. 

In conclusion, both our numerical data and mean field considerations suggest that the 
Z(2) and U(l) Yukawa models in 2d have asymptotically free Yukawa coupling and exhibit 
dynamical fermion mass generation also for arbitrarily strong scalar field selfcoupling, when 
the usual perturbative and 1/Np expansions are not applicable. They are thus quite similar 
to the corresponding Gross-Neveu models. However, at present we cannot say whether 
they belong to the same universality class, i.e. what the precise form of the scaling law is. 
Our precision is also not yet good enough to compare with the results of the mass gap 
calculations ||. Further study of these questions, as well as of the transition from the 
Gross-Neveu-like behaviour to the NLcx or spin model cases is needed. Finally, we have not 
yet studied the properties of the 2d lattice Yukawa models at larger negative k, where the 
antiferromagnetic scalar field coupling competes with the Yukawa coupling which supports 
ferromagnetic order. 
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Figure 1: The numerical data for the fermion mass arriF on finite lattices in the U(l) Yukawa 
model plotted against 1/y 2 , where y is the Yukawa coupling constant. The dotted, dashed and 
full lines are the mean field predictions on the L 2 lattices, L = 16, 32 and 64, respectively. The 
statistical error bars are omitted to make the plot tidy, they do not exceed the symbol size. In this 
figure k = A = 0, which means that the U(l) Yukawa model is identical to the chiral Gross-Neveu 
model. 



Figure 2: As in fig. 1, but now for k = 0.2 and A = 0. Here both the 1/Np expansion and the 
MFA are as well applicable as in the GN case. 



Figure 3: As in fig. 1, but now for k = and A = oo. The 1/Np expansion is not applicable but 
the MFA still describes the data very well. 



Figure 4: As in fig. 1, but now for k = 0.2 and A = oo. The 1/Np expansion is not applicable and 
also the MFA does not describe the data very well. The dotted-dashed line represents a fit to the 
64 2 data (circles) by means of the generalized gap equation ( |4.1| ). 
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